Macroscopic Quantum Mechanics in a Classical Spacetime 
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We apply the many-particle Schrodinger-Newton equation, which describes the co-evolution of an many- 
particle quantum wave function and a classical space-time geometry, to macroscopic mechanical objects. By 
averaging over motions of the objects' internal degrees of freedom, we obtain an effective Schrodinger-Newton 
equation for their centers of mass, which are degrees of freedom that can be monitored and manipulated at 
the quantum mechanical levels by state-of-the-art optoemchanics experiments. For a single macroscopic object 
moving quantum mechanically within a harmonic potential well, we found that its quantum uncertainty evolves 
in a different frequency from its classical eigenfrequency — with a difference that depends on the internal 
structure of the object, and can be observable using current technology. For several objects, the Schrodinger- 
Newton equation predicts semiclassical motions just like Newtonian physics, yet they do not allow quantum 
uncertainty to be transferred from one object to another through gravity. 



Introduction. — Prescribing and testing modifications to non- 
relativistic macroscopic quantum mechanics due to self grav- 
ity has been one, although minor, approach towards the explo- 
ration of nontrivial effects of quantum gravity. Apart from the 
standard formulation of linearized quantum gravity [ 1 ] , which 
seems rather implausible to test in the lab, several types of 
schemes have been proposed [2-12]. The first type is grav- 
ity decoherence [2-6, 13-15], where the existence of grav- 
ity is conjectured to introduce decoherence to macroscopic 
quantum superpositions. The second type modifies canonical 
quantization motivated by the existence of a minimum length 
scale [7-9]. A third type, which will be the subject of this 
paper, is often referred to as semiclassical gravity [10-12]. 
As originally suggested by Moller [10] and Rosenfeld [11], 
spacetime structure might still remain classical even if it is 
sourced by matters which have quantized motion. The idea is 
to impose the following condition {G - c - 1): 



(1) 



where G^,, is the Einstein tensor of the (3-i-l)-D classical 
spacetime, 7),v is an operator for the energy-stress tensor, and 
\4i{t)) is the wave function of all matter in the universe, which 
evolves within this classical spacetime. 

Many arguments exist against semiclassical gravity. Some 
rely on the conviction that a classical system (gravity) can- 
not properly interact with a quantum system (electroweak plus 
strong interactions) without creating contradictions. Others 
are based on "instrinsic" mathematical inconsistencies, the 
most famous one between Eq. (1), state collapse, and the 
Bianchi Identity [16]. Towards the former type of argument, it 
is exactly the aim of this paper to work out the effects of grav- 
ity being classical on the quantum mechanics of macroscopic 
objects. Although we will find them counter intuitive, they do 
not seem to us as completely dismissible right away. In fact, 
we shall find these effects "right on the horizon of testability" 
by current experimental technology. Towards the latter type of 
arguments, we shall remain open minded regarding the possi- 
bility of getting rid of state reduction while at the same time 
avoiding the many-world interpretation [17]. 



We will consider the non-relativistic version of Eq. (1), 
the so-called Schrodinger-Newton (SN) equation, for macro- 
scopic objects each consisting of many particles, and show 
that within certain experimental parameter regimes, the 
center-of-mass (CM) wavefunction approximately satisfies 
the following nonlinear Schrodinger-Newton equation (with 
Ti = 1): 
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Here (x) = (^\x^) is the expectation value of CM position; 
Wc is the eigenfrequency of CM motion in absence of self 
gravity; we shall refer to C as the SN coupling constant, and 
wsN = ^CjM the SN frequency. 

For a single macroscopic object prepared in a squeezed 
Gaussian state, Eq. (2) leads to a different evolutions of ex- 
pectation values and quantum uncertainties, as illustrated in 
Fig. 1. Such a distinctive deviation could be tested by op- 
tomechanical devices in the quantum regime [ 1 8-2 1 ] . For two 
macroscopic objects interacting through gravity, we show fur- 
ther, using the two-body counterpart of Eq. (2), that classical 
gravity cannot be used to transfer quantum information — al- 
though experimental demonstration of this effect will be much 
more difficult than the modification to single-object dynamics. 

Here we emphasize that it is not our aim to explain the col- 
lapse of quantum states using the SN equation, as has been 
attempted in the literature [22-24] . In fact, we will avoid alto- 
gether the regime in which wavefunctions can be highly dis- 
torted [25-27] by the nonlinearity of the SN equation, or other 
proposed quantum-gravity effects, and constrain ourselves to 
Gaussian states whose evolutions only deviate very little from 
predictions of standard quantum mechanics. The experimen- 
tal tests we propose will require less isolation from the envi- 
ronment, for much shorter periods of time, although the price 
we pay is higher level of dependence on the models we are 
testing, and less dramatic symptoms for those models. 

Many-particle SN equation. — For n particles, we denote 
their joint wave function as i^(f, X) with 3«-D vector X = 
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FIG. 1. (Color online). By contrast with the standard Schrodinger 
evolution (left panel), the covariance ellipse of a Gaussian state, un- 
der the Schrodinger-Newton Equation (right panel), rotates at an an- 
gular frequency io^ = (lj^ + w^j^)''^ higher than io^, at which the 
vector ((x), (p)) rotates. 



(xi, ■ ■ ■ , x„) and the 3-D spatial coordinate of ^-th par- 
ticle. From Refs. [10, 11], the co-evolution of the (3-i-l)-D 
spacetime background with coordinate (f , x) and the «-particle 
wave function is given by [cf. Eq. (1)] 



d'"XT.y(x,Xmt,X)\\ 



Q^Yj^Uu+mlm, X) 



(3) 
(4) 



where Da is the (curved-spacetime) D'Alembertian with re- 
spect to the A:-th spacetime coordinate and is the mass 
of the A:-th particle. In the non-relativistic limit, Diosi and 
Penrose [5, 22] factored out a fast-varying phase, \j/(t, X) = 
e-'Zt'«A'^(f^ X), and obtained 

id,ip = Yj [-'^i/(2'«^) + m Uit, x,)/2] ip + ViX)<f , (5) 

k 

where V(X) is the energy for non-gravitational interactions, 
and the Newtonian potential U{t, x^) is given by 

V-f/(f, x) = 4;r^ J d^"X\ip{t,X)\-mj6ix- Xj) . (6) 



Equations (5) and (6) are still not concrete enough for exper- 
imental studies, because we cannot separately probe or drive 
the motion of each particle of a macroscopic object. In op- 
tomechanical devices, a light beam often probes (and hence 
acts back onto) the average displacements of atoms within the 
first few layers of the reflective coating of a mirror-endowed 
mechanical resonator Motion of this effective surface loca- 
tion can often be well-approximated by the CM motion; the 
error of this approximation is referred to as the "internal ther- 
mal noise", and has been shown to be suppressible below the 
quantum level of CM motion [28-30]. 

Separation of scales. — Before actually deriving the CM equa- 
tion of motion, let us survey the time and length scales present 
in our system, which are critical for separating the CM mo- 
tion from the internal motions. To be specific, let us con- 
sider a macroscopic object (crystal) with equal-mass particles 
(atoms, mass m) positioned on a uniform 3-D lattice. 



For time scales, the CM motion is determined externally 
by the measurements we chose to perform. In this paper, 
we consider motions from Hz to kHz scale. By contrast, 
the internal motion — the oscillation of the nuclei around 
their equilibrium positions on the lattice — takes place at the 
Debye frequency aio of the material, which is around tens 
of THz — much faster than the CM motion [31]. The co- 
herence time of CM motion can be made very high, up to 
TCM ~ Q/<^CM ~ lO'^^/wpiyi, where Q is the mechanical qual- 
ity factor and we have used the fact that Q wcm ~ 10^^ ~ 10^"^ 
for low-loss materials [20]; tcm can be minutes or even hours, 
and this is orders of magnitude longer than the coherence time 
Tint of the internal motion, which is usually less than /is, even 
under cryogenic temperature [32]. 

For length scales, the typical CM motion, during a 
quantum-limited measurement process, is characterized by 
its zero-point motion, and Axcm(10"^^-10 kg, 1-1000 Hz) ~ 
10 '^-10 '^m. In comparison, the zero-point motion for an 
internal degree of freedom (DOE), i.e., a nucleus, is given by 
AjCint ~ y]til{2m Uind where m is the mass of the nucleus, usu- 
ally between 10 and 100 times the mass of the proton (~ 28 for 
a Silicon crystal). Since the typical frequency of the internal 
motion Wint is close to the Debye frequency wd ~ 10'^ s"', 
we have IS.Xm ~ 10 '^m, which is much larger than Ajccm- 
Therefore, the quantum mechanical CM motion we impose 
externally is tiny compared with internal motions. 

SN equation for the CM. — Keeping the above separations 
of scales in mind, we derive the SN equation for the CM. 
Suppose we have a crystal with « atoms, then the CM is at 
xcM = (l/n)ZijiXi: and the internal motion of the ^-th parti- 
cle with respect to the CM is defined as = Xk - xcm- In 
standard quantum mechanics, given a Hamiltonian that only 
depends on the separation of particles, the CM and the inter- 
nal DOEs are separable: (f{t, X) - '^cuiU ^cyiY^mtit, Y), with 
3(« - 1)-D vector Y s (yi, ■■■ , y„-i). The two wavefunctions 
evolve independently when the external force acts directly on 
the center of mass [through y(x)]: 

id,^'cM{t, X) = [-V-/(2M) + V{x)\ ^>cm{U X) , (7) 
!5,>Pi„,(f,Y) = //iA(r,Y). (8) 

Here i/int acts only on the the internal DOEs. 

In the case of the SN equation, however, the Newtonian 
potential U in Eq. (5) does not allow a strictly separable evo- 
lution of ^FcM and ^int. Nevertheless, as we will show, the CM 
and the internal DOEs can be approximately separable at time 
scales long compared with the coherence time of the internal 
DOEs, but much shorter than the CM dynamical time scale. 
Suppose at time t, we have a product of pure states of CM and 
the internal DOEs, and let us consider an evolution from t to 
/ H- Af with increment Af satisfying 



TCM » UcM S Af » Tint » 



(9) 



therefore focusing on the evolution time scale of the CM. 
Within Af, the internal DOEs suffer from decoherence. To 
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accommodate this fact, let us evaluate the state evolution by 
using density matrices. At moment f, the joint density ma- 
trix is given by p(f) = \Lp(t)){ip(t)\ = pcm(0 ® Pint(0, where 
we have introduced pcm(0 = I*I'cm(0><*I'cm(OI and pimCO = 
l*I'int(0){*ifint(OI- Sincc wc are mainly interested in the effect 
of the Newtonian potential U, we move into an interaction 
picture in which U is the only interaction Hamiltonian. Up to 
the first order of Af, we have 



pi(t + At) -p,(t) 



At'[u,{t'),pm (10) 



with p/(f) - p{t) - pcm(0 ® Pint(0 initially. Here the New- 
tonian potential operator f//(f) is a simplified notation for 
f//(r, Y(f),xcM(f)), defined as [cf. Eqs. (5) and (6)] 



U:(t) = m 



'2/ 

k.j ^ 



d^x 



Tr{p/(f)5[y/f)+XcM(0-x]} 
|x - y*(f) - xcm(OI 



(11) 



Note that yk(t) and xcM(f) are operators in the interaction pic- 
ture with time evolution under the free Hamiltonian and the 
influence of thermal decoherence. 

Because A? » Tint, the internal DOFs undergo strong deco- 
herence and we can therefore assume the internal motion Y(f) 
in the interaction picture is ergodic, and the time average can 
be approximately by the ensemble average: 

df';7/(f') = Tri„,[p™f//] + ee„-. (12) 



-f 

Af J, 



The ensemble density matrix of the internal motion p™, dif- 
fering from pi{t), is equal to the one averaged over thermal 
equilibrium. The error term ee„ in this approximation is sup- 
pressed by two factors: (i) the small ratio between the coher- 
ence time Tint and Af, which breaks fen into many incoherent 
contributions in time, and (ii) the large number of independent 
atoms, which further breaks e^n into many independent con- 
tributions from different atoms. Appendix B provides more 
details for estimating e^n- 

Under assumption ( 1 2), at the time scale we choose, the CM 
and the internal DOFs will remain separable; the CM state 
will remain pure up to the decoherence time tcm, while the 
internal DOF will undergo decoherence at a much faster time 
scale. After singling out the CM, we obtain the SN equation 
for the CM in the interaction picture. 



/5,'PcM=Tri„,bS™C///2]^CM. 



(13) 



In Appendix A, we show that when Axcm ^ AjCint, keeping 
the leading order in xcm, we can write 

TrintbS™ U,] ^ £"(0)/2 [<4^) - 2(icM)A-CM + 4m] (14) 

after having restricted center-of-mass motion along one par- 
ticular direction. Here (xcm) = (*I'cmI-?cmI*I'cm); &" is the 
double (directional) derivative of the Newtonian potential en- 
ergy £(x) between the mass distribution according to the in- 
ternal motion and the same distribution but translated by x, 

eint(x')£'int(x") 



£(X): 



Ix' - X" -I- xl 



(15) 



with gint(x) the matter density at position x, given by: 

eint(x) = Yj r '^""''^ Tr [<Y|p™ |Y)5(x - y,)] . (16) 

Neglecting the c-number term {x^y^ - (j^cm)^ as it only adds 
an overall phase to the wave function and putting back the free 
part of the Hamiltonian for the CM, we obtain Eq. (2) shown 
at the very beginning with SN coupling constant defined by: 



C = &"(0)/2 . 



(17) 



Estimates for wsn- — Let us now estimate the magnitude of 
a>SN- Suppose we have an object with mass M and size L; 
assuming an absolutely homogeneous mass distribution, we 
obtain 



6^°"" ^ 3GM^/{2L^) ^ 3GMpo/2, w^™ ^ ^/3G^, (18) 

where we put back the Newton's constant G and po is the den- 
sity of the homogenous material. This is a typical order of 
magnitude obtained in the literature. However, according to 
results obtained so far in the paper, we must account for the 
concentration of material near the equilibrium positions of the 
nuclei, with Axint ^ ^lattice, we can roughly model the ob- 
ject's matter distribution according to pint as a square lattice 
of solid balls with radius Axi^t and constant density — with 
lattice spacing equal to aiattice- It is easy to argue that the self 
gravitational energy S of such a lattice is dominated by the 
sum of the self energies of each ball, which gives 



(19) 



In other words, we have an amplification factor A » 1 for 
the SN constant. For Silicon at T = 100 K, we already have 
H(Dd ~ IkgT and the variance of nucleus thermalized motion 
is about the same as Ax^^^. Noting further that aiattice ~ 5 x 
lO^'^m, Po ~ 2.7g/cm-', and Axint ~ 5 x 10"'~m, we have 



A ~ 10" 



crystal 



'0.4 s" 



(20) 



We shall use this (^5^ as our baseline estimate for ^l>sn 
throughout the rest of the paper 

Evolutions of Gaussian States. — As one can easily prove, an 
initial Gaussian state evolves under Eq. (2) will remain Gaus- 
sian. Here we write down the self-contained evolution equa- 
tions for first- and second-moments of x and p, which com- 
pletely determine the evolving Gaussian state: 



Md{x)ldt = d{p}/dt, d{p}/dt = -Mulix) , 



(21) 



V„- = 2y,„/M, Vp 



-2M(cjI + cjI^)V,p , (22) 



(23) 



For covariance we have defined Vab = (AB + BA}/2 - {A}{B}. 
Equation (21) indicates that expectation values of x and p 
evolve the same way as a standard simple harmonic oscilla- 
tor with angular frequency w^; this means any semiclassical 
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measurement of the oscillator on expectation values of x and 
p will confirm classical physics. Evolutions of second mo- 
ments, i.e., the covariance matrix, which represent quantum 
uncertainties, however, are modified; they evolve the same 
way as a harmonic oscillator with a modified frequency: 



+ Ll)\ 



SN ■ 



(24) 



as illustrated in Fig. 1 . In Appendix C, we show that the differ- 
ence between first- and second-moment evolution may show 
up in the output spectrum of an optomechanical system in the 
quantum regime. In order to measure the effect of the SN 
term, one must be able to distinguish between oj^ and Wq, by 
requiring oj^-oJc > ciidQ, where Q is the quality factor of the 
mechanical oscillator. This requires 



Q I (wsn/wc)-. 



(25) 



which is experimentally achievable if wsn ~ 0.4 s uqu ~ 
2n X lO^Hz, and e ; 2 X 10** [33]. 

SN equationfor two macroscopic objects. — If gravity is intro- 
duced as an interaction potential between two objects which 
are confined within potential wells with common eigenfre- 
quency wq, and moving along the same direction as the sepa- 
ration vector L between their equilibrium positions, we have 



V - CiiixcM, - XCM2) 



with 



= 2^' ^'^ 



|Xi - X2 - L| 



~(2)/ 



(26) 



(27) 



where p[ll and p'^l are ensemble-averaged mass distribution 
densities for objects 1 and 2, respectively. This weak coupling 
makes differential mode of the two CM motions oscillate at a 
slightly higher frequency — thereby allowing quantum state 
to slosh between the two objects, at a frequency of A = \(l)+ - 
(^-\ = C\2/{2Mcl)()). In order for information to successfully 
transfer, one requires both objects to be in equilibrium with a 
heat bath with kgT i ftQ^, and A J coq/Q, which requires 



Q > MullCn . 



(28) 



Suppose we instead use the SN equation for two macroscopic 
objects in this scenario. Following the same procedure as 
above, we obtain the SN equation for the CMs of two objects: 



(29) 



Here H\i 22 are the sum of the free Hamiltonian and the SN 
interaction term for each individual CM [as in Eq. (2)] , £12 
is the zero point of the mutual Newtonian interaction energy 
and H\2 is given by, 

Hn = C12 [(xcMi - {xcM.)f + {xcM2 - <^CM, )f\ /2. (30) 

Equation (30) makes sure that only expectation value of posi- 
tion gets transferred between the two objects (same way as in 



classical physics), while quantum uncertainty does not trans- 
fer from one object to the other In retrospect, the expecta- 
tion value gets transferred because the semiclassical limit of 
the SN equation is simply Newtonian physics; quantum un- 
certainty does not get transferred because classical space-time 
is incapable of storing, let along passing on, quantum uncer- 
tainty. 

However, because C12 ~ GM^/L^ i GMpo <sc C, the lack 
of the amplification factor A [Eq. (19)] in C12, makes the Q 
required [Eq. (28)] for testing the inability of classical gravity 
to transfer quantum information is much larger than the the 
threshold for testing the self-gravity effect. We will be forced 
to require ~ 2n x I Hz, if g ~ 2 x 10*^, which seems very 
challenging for ground-based experiments, yet might not be 
impossible [34]. 

Discussions — In this paper, we have taken the liberty to as- 
sume that spacetime might be classical [12]. In our opinion, 
due to the lack of precision experimental tests on the quan- 
tum coherence of dynamical gravity, semiclassical gravity is 
still worth testing. In fact, our calculations show that detecting 
the consequences of classical gravity using quantum mechani- 
cal phenomena requires carefully designed experiments using 
techniques that have not been available until recently. 

A concrete mathematical formulation of non-relativistic 
semiclassical gravity, namely the SN equation, combined with 
basic knowledge of condensed-matter physics, has revealed 
self-gravity effects in macroscopic objects that can be tested 
much more easily than performing the more obvious (but less 
model dependent) experiment of demonstrating impossibility 
of information transfer Thanks to the amplification factor 
A, this "much easier test" seems only moderately challeng- 
ing with current technology. 

We speculate that the rate gravity decoherence should also 
be expedited by the same A'^^ ~ 1000, because gravita- 
tional self energy is also used as an indicator of decoherence 
time [5]. However, due to the lack of a widely-accepted mi- 
croscopic model of gravity decoherence, this only makes it 
more hopeful to make experimental attempts, but would not 
enforce a very powerful bound if decoherence were not to be 
found. 

Since our classical gravity model literally requires the exis- 
tence of a global wave function of the universe that does not 
collapse, in the arguably unlikely case the proposed experi- 
mental test shows any positive result, we will open up the op- 
portunity of investigating the fundamental nature of quantum 
measurement using gravity. 
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APPENDICES 

A. Details in calculating Trimlp™ f/;] 

In this section, we will fill in details that would lead to 
Eq. (2). Let us first define 



= yt +xcM, k=l,...,n. 



(A.l) 



Here are positions of the ^-th particle, xcm is the center- 
of-mass position, while are internal motions, with the con- 
straint that ^11=1 yi< ~ 0- ^'^^ ^ density matrix p that is fac- 
torable into p - p™ ® pcM, the CM and the internal motions 
are independent from each other For this p, we can easily 
obtain the matter-density distribution 



k=i 



k=l ^ 



d"" Yc/-'xcMPmt(Y)pcM(xcM)5^^^[x - (yA + Xcm)] , 

(A.2) 

where for convenience we have defined 

Pint(Y) = <Y|p™ |Y> , pcm(xcm) = (xcmIpcmIxcm) ■ (A.3) 

In a similar way, we can also define matter-density distribu- 
tions given by the internal motion alone (assuming xcm = 0), 
as well as the center-of-mass matter-density distribution: 



N 

eint(x)=^m <f"-^Ypi,A^)6^^\^-y,]. 
I — 1 «^ 



k=l 

PcmCx) =mpcM(x) . 
Together, Eqs. (A.2) and (A.4) gives 



etot(x) 



PintCx - y)gcM(y)d y , 



(A.4) 
(A.5) 

(A.6) 



which is rather intuitive. We can then write 

N N 



etot(x) 



k=l 



This leads to 



|x - y^t - xcmI 



V r £'tot(x)5*-''(z -y*) ,1 
j-^J \x-z- XcmI 



xd^z (A.7) 



r N 

I k=l 



-! 



gi.t(X-y)ge.(y)e.n.(z)^3^^3y^3^ 

|x - z - XcmI 

r ,3 / N rj3 ,3 emt(x)£'int(z) 

■■ I d y£)cm(y) I d xd z- 



\x-z + y- XcmI 



(A.8) 



At this stage, if we define 

&(x) = J d'yd'z 



3.,J^ eint(y)eint(z) 

|y-z + x| 



(A.9) 



this is the mutual gravitational potential energy between a 
matter distribution Qi^i and an identical one translated by x. 
If translation is along one direction, with a magnitude much 
less than the scale in which £ changes, then 



fi(x) * £"(0)xV2 + £(0) 



(A. 10) 



The constant energy piece does not affect wave function evo- 
lution except adding an overall phase — which can also be 
absorbed into the definition of the internal wavefunction. This 



means 



t^mxlPZUi'] = \£"{0) j (y-xcM)Vm(y)t/V 



ifi"(0) [<4m) - 2<icM):fcM + 4m] ■ (A.l 1) 



B. Estimation of the Fluctuation Term 

In this section, we estimate the error of the fluctuation 
which appears in Eq. (12); 



At J, 



dt'Uiit') =u + ee„. 



with 



U =trir 



pZt,'nU(p,x,) 



k=l 
-/+A/ 



Gnrn^ . 2 



Axr 



3 "-^CM 



^Xt dt'Y,[Uk(t')-{Uk(t')}] 

^' k=i 



(B.l) 

(B.2) 
(B.3) 



Here the (...) in eg,, represents ensemble average. It would 
not be difficult to estimate, by squaring Eq. (B.3) and taking 
expectation value, that 



{elu) ~ m J —(AC/) ~ 

V At Axint V nAt 



(B.4) 



Here t» is the coherence time of Uk{t), AU is the likely range 
in which U can fluctuate, while the number of particles « en- 
ters this way because motions of atoms are largely indepen- 
dent from each other This means we have 



IAxcm/ V nAf 



(B.5) 



Here the separation of time scales, plus the large number of 
n, can often overcome Axint/AxcM factor and make e^j 
u. This means we can ignore the error introduced by using 
ensemble average. 
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C. Effective Heisenberg Equations of Motion and Coupling witli 
Optical Field 

The fact that Gaussian states leads to Gaussian states en- 
courages us to look for effective Heisenberg equations of mo- 
tion, which will at least be valid for Gaussian states. It is easy 
to find that 

dx/dt^p/M (C.l) 
dp/dt - - MloIx + MljI^^{x - {x}) (C.2) 

will give the same set of first- and second-moment equations 
of motion. Note that in the Heisenberg picture, the initial state 
of the oscillator remains constant. 

Let us now consider a more realistic scenario, in which the 
oscillator is damped with decay rate driven with classical 
thermal noise and other classical driving; we also consider 
using light to sense the position of the mirror, in which we also 
suffer from sensing noise. The entire process can be described 
by the following set of equations: 

dx/dt^p/M, (C.3) 
dp/dt - - M[(L)lx + 2y„ip + colf^ix - {x})] + aa + f , (C.4) 
bi —02 + n + ax , bi — ai . (C.5) 

Here y,,, is the rate of damping, a is the optomechanical cou- 
pling constant, / is classical driving (we ignore the tiny quan- 
tum contribution from the thermal force). We have used 01,2 
to represent the quadratures of the incoming optical field, and 
^1,2 those of the out-going field. We have used n to denote 
classical or quantum sensing noise. 

This set of equations can be solved first for (i) by taking 
the expectation value of the first two equations, and then insert 
this back to obtain the entire solution. If we define 



1 



(C.6) 



M{co^ + 2iy„,cu - tol) 

" 'l^TTTT^ ^ ' ^1 " a/'^c (C.7) 

M(a)^ + 2iy„,(jL) - Wq) ^ 



then we can write, in the frequency domain, 
^2 = ^2 + aXgh +Xof + n . 



(C.8) 



Because both Xg and xo in the time domain are Green func- 
tions of stable systems, Eq. (C.8) represent the steady-state 
solution for the out-going field, which is only determined by 
the in-going optical field and the classical driving field — ini- 
tial states of the mechanical oscillator does not matter [similar 
to the caseofRef. [35, 36]]. 

Equation (C.8) carries the separation between classical and 
quantum rotation frequencies in the previous section (Fig. 1) 
into the frequency spectrum of our measuring device: quan- 
tum radiation-pressure noise spectrum in the output port of 
the continuous measuring device is the same as an oscillator 
with frequency Wq, and therefore peaks around Wq — while 



classical noise follows that of an oscillator with frequency Wc, 
and peaks at oj^. In order to look for such a signature, we will 
need classical force noise to be comparable in level to quan- 
tum noise, and have the two peaks to be resolvable, 

Wq - J Jm (C.9) 

which means 

Ql(a>,la),nf (C.IO) 
where Q is the quality factor of the mechanical oscillator 

D. SN Equation for Two Macroscopic Objects 

Following the analysis in Appendix. A, here we deal with 
two macroscopic objects. For simplicity, let us assume the 
two objects as identical, and define 

4''=yf'+4M' k^\,...,n, 7=1,2. (C.ll) 

As the distance between the two objects is macroscopic, their 
wavefunctions have no overlap. We also have 

eSCx) = J eL'iCx - y) d^Wy , / = l, 2 . (C. 12) 

The total SN term in the joint SN equation will be 

2 n ^ (I) 



J— 1 1 f f_ 1 ; — 1 ^ A 



(x) , 

d^x (C.l 3) 



/=1 <:=1 I.J=\ k=\ ^ I* '1 

It is easy to identify contributions to the self-SN terms: 



iiJ |x-xl'>| £f J |x- 



^>(2)| 



(C.14) 



The mutual SN terms can be evaluated via ensemble averaging 
over p^" - P- , ®p, : 

int int int 



tr;, 



iiJ |x-4 



~£l2(0) + £'12/2 [-^CM, - 2xcMi<-XCM,> + <^CM,)] (C.15) 
~£l2(0) + £'12/2 [-^CM, - 2XCM2(-^CM,> + <^CM|)] (C.16) 

where fin is given by 



fii2(x)= dyd^z 



,3 ,3 £'!::(y)eL'V) 



|y - z + L + x| 



(C.17) 



